Controlling spin motion and interactions in a one-dimensional Bose gas 
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Experiments on ultracold gases offer unparalleled opportunities to explore quantum many- 
body physics, with excellent control over key parameters including temperature, density, 
interactions and even dimensionality. In some systems, atomic interactions can be adjusted 
by means of magnetic Feshbach resonances, which have played a crucial role in realizing 
new many-body phenomena. However, suitable Feshbach resonances are not always avail- 
able, and they offer limited freedom since the magnetic field strength is the only control 
parameter. Here we show a new way to tune interactions in one-dimensional quantum 
gases using state-dependent dressed potentials, enabling control over non-equilibrium spin 
motion in a two-component gas of 87 Rb. The accessible range includes the point of spin- 
independent interactions where exact quantum many-body solutions are available and the 
point where spin motion is frozen. This versatility opens a new route to experiments on 
spin waves, spin- "charge" separation and the relation between superfluidity and magnetism 
in low-dimensional quantum gases. 



Advances in optical and magnetic trapping of ultracold gases have played an essential role 
in opening up novel avenues in quantum many-body physics by providing experimental access 
to new physical regimes In particular, one-dimensional (ID) quantum gases, created using 



optical lattices or atom chips, exhibit a surprisingly rich variety of regimes not present in 2D or 
3D JjJ]. For example, a ID Bose gas becomes more strongly interacting as the density decreases. 
Furthermore, the many-body eigenstates and thermodyna mic p roperties of these ID systems can 



often be described using exact Bethe Ansatz methods 



theory and experiment are possible 



9. 



1CH14I]. and direct comparisons between 



15H17I]. Adding the possibility to d yna mically control 



16|, 



20] , there is now a 



the strength of atomic interactions, for example via Feshbach resonances 
strong impetus to extend these experimental and theoretical studies to non-equilibrium dynamics. 
Spinor quantum gases offer the opportunity to study the interplay between internal (spin) 



and external (motion) degrees of freedom 



161, 



13, 



21 



271 ] . In this context, strong candidates for 



experiments are the two magnetically trappable clock states in °'Rb 



23 



251 ]. in part 



because 



they experience equal trapping potentials and have nearly spin-independent interactions 



2a-l30(]. 



The drawback is that no convenient Feshbach resonances are available for these states, preventing 



2 



precise control of the three relevant (inter and intranstate) interaction strengths. 

In the two-component ("spin- 1/2") ID Bose gas, the presence of spin-independent (symmetric) 
interactions is of particular interest. For all interaction strengths (weak and strong) the dispersion 



23 



311 ] . and the low-energy spin velocity vanishes. As a 



relation of spin waves is quadratic here 

consequence the usual Luttinger- liquid description [s , If], 32], 3s| c 
precisely the point where exact Bethe Ansatz methods can be used Q, Furthermore, it is the 
point where buoyancy effects vanish and in the weakly interacting (mean-field) regime it also lies 



3l|. 



on the border that separates miscible and immiscible regimes of binary superfluids 

We show that radio-frequency-dressed potentials on atom chips offer a new way to tune the 
effective interactions in ID and to control spin motion. We make use of the fact that, for elliptical 
rf polarizations, different hyper fine states experience different dressed potentials, allowing for state- 



dependent manipulation [34(. Here we exploit the dependence of the ID coupling strength on the 
transverse confinement frequency co± Q]. State-dependent optical lattice potentials have previously 



found use for spin-dependent transport and entanglement of atoms 



35 



3d ]. More recently, state- 



dependent microwave dressing was used to generate spin squeezing in 3D Bose-Einstein condensates 



37 



and state- 



by varying the wavefunction overlap for two hyperfine states to control collisions 
dependent potentials created by combining an optical trap with a magnetic field gradient were used 
to obtain record low spin temperatures via spin gradient demagnetization cooling of a quantum 



gas 



391 ]. By tuning the transverse confinement for the two states independently through the rf 



polarization and amplitude, we show that it is possible to control the interactions in a state- and 
time-dependent manner. Suddenly changing interactions, combined with the state dependence of 
the axial trapping then results in dynamical evolution in the spin degree of freedom. In particular, 
we are able to tune to (i) the point where the spin motion is frozen, and (ii) the point where the 
ID interactions become spin- independent. 

We first discuss our results on the one-dimensional non-equilibrium dynamics for state- 
independent potentials, highlighting the importance of small differences in interaction param- 
eters. The starting point of our experiments is a nearly-pure ID quasi-condensate in the 
|1) = \F = l,mf = —1) state of 87 Rb in a highly elongated magnetic trap created by an atom 
chip (see Methods). From this initial state, we induce a sudden transition to a coherent super- 
position of the |1) and |2) = \F = 2, mt = 1) hyperfine states via a two-photon pulse, effectively 
creating a spin- 1/2 system |2£ 



\ Q] • The resulting non-equilibrium situation is allowed to evolve for 
a variable hold time. Subsequently, we directly image the longitudinal distributions, and obtain 
the linear densities n\ and ri2 of the two states along the length of the trap. 
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FIG. 1: [Color] One-dimensional spin dynamics and total density after a sudden transfer of internal-state 
population, for the case of state-independent trapping potentials. Shown are spin polarization (rii — ri2, 
left) and total linear density (n^ + n 2 , right), as a function of axial position and time after the transfer. 
Top: experiments, bottom: corresponding simulations resulting from integration of two coupled ID Gross- 
Pitacvskii equations (GPE). The spin polarization data clearly shows how n<i is focused towards the center 
(blue), while n\ moves towards the sides (red); the total density shows little dynamics. Differences between 
experiment and simulation can be explained by the limited optical resolution of our imaging system and a 
small tilt of the trap, leading to a slight spatial asymmetry in the experiments. 



In figure Q] we present measurements of the evolution of spin polarization (ni — 71-2) and the 
total linear density (n\ + n.2) as a function of hold time. The spin pattern shows clear dynamical 
evolution [fig. [U^a)] whereas the total density remains approximately constant with no significant 
dynamics [fig. QJb)]. The spin dynamics can be interpreted as a "focusing" of state |2) in the 
presence of state |1), resulting in a negative spin polarization (ri2 > rii) toward the center of the 
trap. 

We find good agreement with the experimental data using the coupled ID Gross-Pitaevskii equa- 
tions (1D-GPE) with solutions also shown in fig. QJc,d). The 1D-GPE is obtained by integrating 



the full 3D-GPE over the transverse ground-state wavefunctions 



411 ]. with interaction parameters 



4 



29)]: a n = 100.4 • a , 



derived from the intra- and interstate scattering lengths taken from ref. 
a22 = 95.00 • ao and a\2 = 97.66 • ao, where ao is the Bohr radius. Generalizing for state-dependent 
harmonic confinement (as will be relevant below) we obtain for the ID interaction parameters 

nn = 2?ia;_|_ ) iaii, 

U12 = AH- j r-Ol2, (1) 

+a;_L i2 ) 

with lo±j the transverse trap frequency for state Similarly we use values for the scaled rate 



29|]. The 



constants for inelastic two-body and three-body losses derived from the 3D values in ref. 
1D-GPE simulations reproduce the features of the experiment, i.e. absence of dynamics in the 
total density and the overall structure of the spin dynamics including the time of maximum state 
separation around t ~ 75 ms. The decay in atom number on a > 100 ms timescale is dominated 
by two-body losses in intrastate interactions and between |2) atoms (712 and 722) 29]. 

The rate of spin focusing/defocusing is critically dependent on the precise differences in ID 
interaction strengths for the respective internal states, a fact that is readily confirmed by changing 
these differences in the simulations. The observed general behavior can be understood as follows: 
in the initial state (an interacting trapped quantum gas in a single internal state in equilibrium) 
the repulsive interactions balance the external confining potential. Suddenly transferring a fraction 
of the population to a second internal state with weaker intra- and interstate interactions results 
in a net contracting force (a confining effective curvature, c% > in equations ([3]) below) on the 
population in this second state that dominates the dynamics in the spin polarization. Because the 
spin-dependent part of the interactions is relatively small, the dynamics in the total density are 
dominated by the (relatively large) average scattering length which remains nearly constant. Hence 
the total density shows only weak dynamics; the focusing in 712 is accommodated by "pushing" n\ 
to the sides (red in fig. [I]). 

We now describe the state-dependent radio-frequency-dressed potentials that we use to control 
the spin motion. We consider near-resonant coupling (huj r f < gpfj,B\B\) of the rf field with tuneable 
polarization determined by the relative phase of two independently controlled rf-fields. A cross 
section of the wire geometry used is shown in fig. [2]^a). The fields originate from direct digital 
synthesis (DDS) supplied currents in two wires neighboring the Z-shaped trapping wire 42|. With 
these two fields we can readily control the ellipticity of the total rf field at the trap position by 
controlling the relative phase (f> of the rf currents in the two wires. This includes linear (horizontal 
and vertical) and circular (cr^) polarizations. 
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FIG. 2: [Color] State-dependent potentials, (a) wire geometry used for state-dependent radio-frequency- 
dressed traps. The static quadrupole magnetic field and the two rf fields are indicated by green, orange 
and purple arrows respectively. The direction of the bias loffe field which defines the quantization axis is 
into the plane of the figure, (b) Trap bottom determined via dressed state rf spectroscopy as a function of 
(ft. Data points correspond to the measured trap bottom for state |1) (red) and state |2) (blue). Solid and 
dashed curves are fits to the data. The dash-dotted green line indicates the fitted trap Larmor frequency 
lol/2-k = 2.25 MHz. 

The corresponding dressed-state potential for state |j) (with j = 1, 2) has the form Vj(x, y, z) = 
((V (x,y,z) - huj rf ) 2 + ^O 2 ) 1 / 2 where V (x, y, z) is the bare magnetic (harmonic) potential. The 



43|, 



44( , which 



state-dependent part of the potential enters through the coupling Rabi frequency re- 
acts to weaken the overall confinement near the trap bottom by an amount given by the dressing 
parameter 5j . Taking the second derivative of the potential Vj around the origin yields new trap 
frequencies, 



5jUJ 2 ± || , where Sj = A/Jfij + A 2 



(2) 



with detuning A = ojl — oj r f and Larmor frequency w^. 

The state-dependent rf p otential is characterized using dressed-state rf spectroscopy with a 
weak additional rf probe 45|, [46|]. Figure [2](b) shows the measured trap bottom as a function of the 
dressing phase (ft, for u r f = 2tt x 2.20 MHz. For (ft = 0.31-rr and (ft = 1.687T the potential is maximally 
state-dependent, corresponding to the pure circular polarizations a~ and cr + respectively (dressing 
only state |1) and only state |2), respectively). The potentials are state-independent for linear 
polarization at (ft = 0,7r (equal dressing of state |1) and |2)). The deviation from a simple sin 2 (</>) 
behavior is due to the wire geometry, as the two rf fields are not quite orthogonal at the trap 
position. A fit to the data (solid and dashed curves) taking into account the wire geometry is used 
to precisely calibrate all parameters of the rf field coupling (see Methods). 



To control the spin motion we turn on the state-dependent dressing, by ramping up the rf 
currents in 2 ms with uj r f = 2tt x 2.20 MHz, directly after preparing the equal superposition of 
|1) and 1 2). The ramp time is slow compared to the inverse Larmor frequency and the inverse 
radial trap frequency, but sudden with respect to any axial motion. We use the two circular rf 
polarizations and various rf amplitudes, corresponding to 0.8 < 5± < 1,^2 = 1 and 0.8 < 82 < 
l,6i = 1. For each time step we extract the widths of the axial distributions in both states. 
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FIG. 3: [Color] Overview of the possibilities of the state-dependent potentials, as a function of the dressing 
parameters (left: varying S2, with Si = 1; right: varying Si with S2 = 1). (a): ID interaction strengths, Uij 
normalised by the bare transverse trap frequency wj_ o- (b) Widths of the distribution at t = 44 ms and 
(c) scaled effective curvature Cj/cq at t = 0. Red indicates state |1) (and un) and blue state |2) (and U22) 
and in (a) U12 is indicated in green. The widths in (b) are obtained by a fit to the experimentally measured 
density profiles (dots) and to GPE simulation (shaded regions). The shaded areas in (b) represent the effect 
of shot-to-shot atom number fluctuations in the experiment. 



Results for the full range of dressing parameters are depicted in figure [3l Figure [3^a) shows 
the calculated interaction strengths taken from equation Q as a function of Si and 82- We have 
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compared the measured widths of the distributions as a function of time with solutions of the 
coupled 1D-GPE. These widths and the corresponding simulations for one fixed hold time of 44 ms 
are shown in figure E^b). The measured widths follow the 1D-GPE simulations closely (taking 
into account the finite optical resolution), with the biggest uncertainties originating from atom 
number fluctuations which cause the peak linear density to vary between 70//m -1 and 100 /xm" 1 
throughout the entire data set (systematic uncertainty shown by shaded regions). The solid vertical 
line at 5± = 0.895 indicates the point where the difference in interaction strengths is minimized 
[hg- EI a )] with un, U22 and u\2 differing by less than 0.05% (100 times reduction in differences 
when compared to the unmodified interactions). These conditions are of interest for comparing to 
Bethe Ansatz solutions which require spin-independent interactions Q, [3] . 

To explain the data we have to consider both the effect of rf-dressing on the collisional interaction 
strengths as well as the state-dependent modification to the axial potential. A simple analytical 
description can be obtained using a Thomas-Fermi description near the cloud center where the 
cloud shape is an inverted parabola. The combination of the state-dependence of the axial trapping 
frequency and of the interactions can then be expressed as a net harmonic potential characterised 
by an effective state-dependent curvature Cj. We solve for the effective curvatures (Fig. [3b) for 
t > in our experiments in terms of 5j, and find 



Here the first term on the right-hand side reflects the modification to the external axial potential 
and the second and third terms deal with the modified interactions Uij . The axial curvature of the 
bare potential is Co = and j3 corresponds to the fraction of the population transferred to 

state 1 2) (/3 ~ 1/2 for our experiments). 

The dashed line at 62 = 0.96 in fig. indicates the point where the difference in interactions is 
compensated by the state-dependent longitudinal potential and cy =02- This point is characterized 
by small equal curvatures of the effective potentials (including interaction energy) for both states 
[fig. Efc)], which result in frozen spin dynamics. These conditions are important for applications 
with on-chip atomic clocks, to minimize inhomogeneous broadening due to mean field shifts. For 
62 < 0.96 the difference in interaction strengths is further enhanced and the time evolution of the 
spin dynamics becomes inverted, with focusing of state |1) while state |2) is pushed outward, as is 
visible in fig. E^b). 

Figure U] shows the full time evolution of the spin polarization for two selected rf-dressing 




(3) 
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parameters. The selected cases are: dressing of state |1) alone ((f) = 0.3l7r, 5i =0.895) [fig. EJa)] and 
dressing state |2) alone ((f) = 1.687T, 52 = 0.96) [fig. IH[b)], corresponding to the intersection points 
in figure 02a) and (c), respectively. Qualitatively state |2) focuses faster with rf dressing applied 
to state |1), when compared to the case of state independent potentials [fig. [Ha)]. Generally, the 
simulated density profiles reveal a rich and dynamic nonlinear evolution of the spin polarization, 
reminiscent of filament propagation in optical systems with competing nonlinearities 4^]. This 
is clearly visible in figure [U^c) for example. This detailed structure depends sensitively on the 
precise values of the dressing. The development and propagation of this fine structure in the spin 
polarization is partially observed in the experimental data, but is not fully resolved due to the 
finite imaging resolution. Convolving the simulated profiles with the point-spread function of our 
imaging system yields excellent agreement with all of the data. With weak dressing of state |2) 
(#2 = 0.96) it is possible to freeze spin dynamics altogether such that the two states maintain their 
overlap and the widths remain constant (apart from a small in-phase quadrupole oscillation and 
decay from state |2)), see Figure [UJb,d). A more quantitative representation of the data, showing 
excellent agreement between experiment and simulation, is given in figure El where the widths of 
the two states are shown for different evolution times. Clearly, the focus point can be identified in 
figure [5]^a) around t = 75 ms and in figure E^b) around t = 30 ms, whereas no focussing is present 
in figure Etc) . 

We have shown that by introducing a small state-dependence to the radial trapping potential 
using rf dressing we can precisely tune the ID interaction parameters in a two-component quantum 
gas by more than 10%, over an experimentally significant range. In our experiments this modi- 
fication competes with the state dependence of the axial trapping and provides a new "knob" to 
control spin motion, leading to tuneable nonlinear behavior. 

Our method can be naturally extended in several ways. For instance, control over the interac- 
tions without the accompanying state-dependence of the axial trapping can be obtained by using 



one-dimensional box-shaped potentials 



421 ]. By introducing an additional displacement of the 



transverse potential in a state-dependent way it is possible to further reduce M12, allowing all three 
interaction parameters to be tuned independently, something that is not generally possible with a 
magnetically controlled Feshbach resonance. 

The observed spin dynamics depend critically on the precise differences in interaction strengths. 
For 87 Rb, the three relevant scattering lengths are nearly equal and therefore weak dressing is 
sufficient to tune the system parameters to the point of symmetric interactions or to where the spin 
dynamics become frozen. Since the rf parameters can be precisely known, such experiments could 
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FIG. 4: [Color] Spatiotemporal behavior of the spin polarization (m — 71-2) following the sudden transfer 
to the state-dependent potentials. States |1) and |2) are indicated red and blue, respectively, (a) shows 
the evolution with rf parameters corresponding to equal inter-atomic interactions (Si = 0.895, ^2 = 1) and 
(b) equal effective potentials (<5i = 1,62 = 0.96). (c) and (d) show the results of 1D-GPE simulations 
corresponding to (a) and (b), respectively. 



also allow precision determination of the scattering length differences. More generally, tuning the 
system parameters around the point of spin-independent interactions strongly affects the dispersion 
relation of the spin excitations 



311 ] . In particular this allows the spin velocity to be tuned around 
zero, providing a new handle for the study of spin waves in one-dimensional atomic gases. 

Tunable interactions in two-component quan tum gases have important applications in the areas 



of spin-squeezing and quantum metrology 



38 



481 ]. and the ability to control spin motion opens 



new avenues for future studies of quantum coherence in interacting quantum systems 
30 



16. 



29. 



I40I |2]. 

at low temperature, and we find that a description based on two coupled ID Gross-Pitaevskii 
equations is sufficient to describe our data. The methods presented here to tune interactions are 
not limited to this regime, however. In particular, we plan to apply these methods to systems with 
stronger interactions (e.g., by lowering the ID density) and with higher temperatures. This will 
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FIG. 5: [Color] Widths of the atomic distributions as function of hold time, for the parameters indicated by 
the three vertical lines in figure [3l Dots are fits to experimental data, lines are results of ID GPE solutions. 
States |1) and \2) are indicated red and blue, respectively. In (a) no rf dressing is applied 5i,2 = 1, as 
in figure [T] (b) and (c) correspond to the data in figure |4] (b) shows the evolution with rf parameters 
corresponding to equal inter-atomic interactions {Si = 0.895,(52 = 1), and (c) equal effective potentials 
{Si = l,S 2 = 0.96). 

provide experimental tests of (and challenges to) more sophisticated theoretical methods, for both 
equilibrium and non-equilibrium phenomena. For instance, it will be possible to experimentally 
explore predictions of the thermodynamic Bethe Ansatz for the two-component Bose gas [bjj] and 



11 



to explore quantum quenches in strongly interacting ID systems by dynamical control over the 
spin-dependent interactions. Finally, we expect that the experimental control over spin motion 
and interactions, as demonstrated here, will benefit the realization of spin- "charge" separation in 



a Bose gas 



23. 



25]. 
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Methods 



Initial-state preparation and coherent spin mixing 

87 Rb atoms in state |1)= \F=l,mf =—l) are evaporatively cooled to quantum degeneracy in a 
highly elongated Ioffe-Pritchard microtrap with trap frequencies of wj_/27r = 1.9 kHz and u\\/2ir = 
26 Hz. The peak linear atomic density is n\ < 100 fim. . In this system, both the temperature 
and chemical potential are small compared to the radial excitation energy (/i, ksT < hu}j_) and the 
dynamics are restricted to the axial dimension (ID regime). A coherent superposition of the |1) 
and |2)= \F = 2^rrif = \) hyperfine states is prepared using a resonant two-photon rf and microwave 
(mw) coupling |2] . The microwave frequency is introduced via an external antenna while the 
rf-field is applied directly to the atom chip wires. The measured two-photon Rabi frequency is 
£l c /2ir = 1.14 kHz, corresponding to a 7r/2-pulse duration of 0.22 ms. This is fast compared to the 
timescale for axial dynamics, but sufficiently slow to prevent radial excitations. Coherence times 
in excess of 1 second have been measured in this setup via Ramsey spectroscopy of dilute thermal 
clouds. 



State-dependent imaging 

The time evolution of the spin distribution is measured by varying the hold time after the 7r/2- 
pulse and performing sequential state-dependent absorption imaging. The atom cloud is released 
from the trap for 1.0 ms time-of- flight to improve the detection efficiency while preserving the 
longitudinal distribution. Atoms in state |2) are imaged directly using absorption on the F = 
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2, F' = 3 transition with an exposure time of 30 /Lis and an optical resolution of 4.2 ftm. State |1) 
is measured by first removing |2) atoms with a resonant light pulse followed by a 1 ms repumping 
pulse from F = 1 to F = 2. The remaining atoms are then imaged in the same way as for state 
1 2). Due to the extra repumping step we find a 20% lower detection efficiency for state |1) and a 
poorer resolution of ~ 8 /xm due to photon recoil, visible in fig. The resulting absorption images 
are integrated along the radial direction to obtain the linear densities m and ni of the two states. 

Characterizing state-dependent potentials 

The potential energy at the trap bottom is characterized by the onset of loss as a function of 
probe frequency which we fit to extract Vj(0, 0, 0). The measured trap bottom varies with rf phase 
between 2250 kHz and 2700 kHz corresponding to a maximum Rabi frequency of fL-/2-7r = 450 kHz. 
A fit to the data (solid lines in figure EJb)) taking into account the wire geometry results in 
an accurate calibration of the key experimental parameters, in particular the Larmor frequency 
Ml = 2.25 MHz, rf field amplitudes b\ = 62 = 0.53 G from the two rf wires and the trap-surface 
distance of 80 /im. 
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